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Abstract. 

In the increasingly active field of optofluidics, a series of experiments involving 
near-critical two-fluid interfaces have shown a number of interesting non-linear effects. 
We here offer, for the first time to our knowledge, an explanation for one such feature, 
observed in experiments by Casner and Delville [Phys. Rev. Lett. 90, 144503 (2003)], 
namely the sudden formation of "shoulder" -like shapes in a laser-induced deformation 
of the liquid-liquid interface at high laser power. Two candidate explanations are the 
following: firstly, that the shape can be explained by balancing forces of buoyancy, 
laser pull and surface tension only, and that the observed change of deformation shape 
is the sudden jump from one solution of the strongly nonlinear governing differential 
equation to another. Secondly, it might be that the nontrivial shape observed could be 
the result of temperature gradients due to local absorptive heating of the liquid. We 
report that a systematic search for solutions of the governing equation in the first case 
yields no trace of solutions containing such features. By contrast, an investigation of 
the second option shows that the narrow shape of the tip of the deformation can be 
explained by a slight heating of the liquids. The local heating amounts to a few kelvins, 
with the parameters given, although there are uncertainties here. Our investigations 
suggest that local temperature variations are the crucial element behind the instability 
and the shoulder-like deformation. 
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1. Introduction 

Optical manipulation of fluid interfaces by means of lasers is a research field in rapid 
growth. Microfluidic applications are already diverse; general accounts can be found, for 
instance, in [1] and j2]. An important advantage of the method is that it is contactless 
and nondestructive, and easily reconfigurable [3H5]. The development of this research 
field has taken place during several decades. Let us briefly mention three milestones of 
this development: 

Our first example is the classic radiation experiment of Ashkin and Dziedzic j6] 

- cf. also Ashkin's extensive reprint volume [7]. Focused light was sent from above 
towards an air-water surface and an outward bulge of the surface of the order of 1 /im 
was observed. The light source was a pulsed frequency doubled Nd:YAG laser, with 
pulse duration 60 ns, peak power 3 kW, and beam waist 4.5 /im. The reason for the 
smallness of the surface elevation was the large air-water surface tension as well as 
the large difference in density between air and water. Theoretical treatments of this 
experiment can be found in [8] and [9] . The second example is the striking experiment 
of Zhang and Chang [10], measuring the deformation of a micro water droplet when 
illuminated by a laser pulse. Typical pulse energies were 100 mJ. Theoretical papers 
can also be found of this effect; cf. |TT] and [T2|[T3]. 

Third — and that is the situation to which we will focus attention in the following 

- is to decrease the fluid-fluid surface tension dramatically by working with a two-fluid 
system of surfactant-coated nanodroplets in oil emulsions near the critical point. A series 
of experimental and theoretical papers have been published by Delville et al. |4~| [T4T[25] 
and by others [2"6||27]. In the vicinity of the critical point the surface tension can be 
made about 10 6 times smaller than the usual air-water tension and the force of gravity 
plays a smaller role since the difference in fluid density is small. The displacement of 
the interface can accordingly be very large, about 70 /xm. 

In the linear regime, when laser power is low, the deformation has been satisfactorily 
explained theoretically using classical electrodynamics — cf., e.g., |26j and [27]. The 
effect is clearly illustrated, e.g. in Fig. 2a of [16]. The phenomenon is sketched in 
Fig. [T^i. However, when the laser power P is increased, typically in excess of 600 mW, 
there occurs a sudden transition into a form illustrated in Fig. b: there is produced an 
elongated lower channel, which we shall refer to as a protuberance. The lower channel 
is narrower than the upper, so that there is a small area of rapidly changing radius, 
a "shoulder", in the displacement. No theory exists to our knowledge for this kind of 
protuberance formation. The effect is definitely of interest to understand, in connection 
with the technique of manipulating soft matter interfaces non-invasively with radiation. 

There are two natural possibilities to explain the instability mechanism behind this 
transition: 

(i) The first possibility is that the reason is of a mathematical nature: the 
governing equation for the surface displacement (Eq. ( |22l) below) may contain the 
protuberance form as one of its solutions, even if physical parameters are assumed 
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(b) Nonlinear deformation with "shoulder" 



Figure 1. Sketch of observed surface deformation in the linear regime, and for higher 
laser powers where the protuberance is observed [16]. 

to remain undisturbed by the laser beam. 

(ii) The second possibility is of a more physical nature: the effect may be due to a local 
change of physical parameters a, p and n caused by the high laser intensity on, and 
in the vicinity of, the central laser beam axis. 

In the following sections we will comment upon the first of these options, and thereafter 
analyze the second one in some detail. First, we will delineate some essentials of the 
electromagnetic theory needed in the problem. 

2. Basic theory 

Assume that laser light comes in vertically from below through medium 2, and becomes 
transmitted into the upper medium 1. The upper medium is the optically denser one, 
so that n\ > n%. We take n\ and n 2 to be real at the actual laser frequency. As for 
gravity the situation is reversed so that the lower medium is the heavier one, p 2 > pi- 
The differences between material constants in the two media are in the present case 
small, of the order of 1%. For the sake of clarity we take the differences to be positive 
quantities, 

An = n 1 -n 2 , Ap = p 2 - p x . (1) 
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Figure 2. The geometry considered: a laser impinges on a fluid-fluid interface from 
below, causing a downward bulge to appear. 



The angle of incidence for the incoming wave with wave vector kj is 9i, the angle of 
transmission (wave vector k t ) is 6 t , and the wave vector for the reflected ray is k r . The 
unit normal vector n is taken to point from medium 1 to 20. For numerical purposes 
later, we shall use numerical values for physical quantities as they appear in |16j . 

Let us delineate how An, Ap and the surface tension a vary with temperature T 
in the vicinity of the critical temperature Tc- First, according to scaling laws 

Ap=(Ap)J^-^-Y , (2) 



. Tc 

3 



where (Ap)* = 285kg/m , with critical scaling exponent = 0.325, and Tc = 308.15 K 
(The separation into two components of the fluid mixture occurs for T > Tc-). 
According to the Clausius-Mossotti relation, An oc Ap, so 



T-T c x p 



An = (An). I J , (3) 
as well, with (An)* = 0.0321. As for the surface tension one has analogously 

with cr* = 1.04 x 10" 4 N/m and v = 0.63 (more details can be found in Ref. [3]). 

X It should be mentioned that these definitions switch I and 2 as compared with earlier works; 

cf. HSU2B1I27I. 
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The pressure difference across the interface z 

1 1 



h due to surface tension is 



p 2 - pi = o — + — 



where R\ and R 2 are the principal radii of curvature. With azimuthal symmetry, as we 
shall assume, 



V% - Pi = - 



ad_ 

r dr 



rh r 



(6) 

0. According to Fig. [2j h < 



with h r = dh/ 'dr. The undisturbed surface is at h = 
for the surface dip. 

Consider now the electromagnetic surface density in the fluid (cf., for instance [9] 
or [28]), 



£ 2 p 



+ 



- 1 d 



£ x H) 



(7) 



We write the constitutive relations as T> = e e£, B = n$H, implying that e is a non- 
dimensional quantity and that all media are assumed non- magnetic. Calligraphic font for 
field quantities indicate real quantities (as opposed to complex field components which 
we will employ later). The first term on the right hand side can be called the Abraham- 
Minkowski term, as it is common for the Abraham and Minkowski energy-momentum 
tensors [9j. The second term is the electrostriction term, important in some experiments 
when the velocity of sound in the fluid is of importance, but not in the present case 
where there is time enough for an elastic pressure to build up to compensate for the 
electromagnetic force [211120] ■ The last term is the so-called Abraham term, also that 
without importance in the present case since this term averages out over an optical 
period. 

What is left is the Abraham-Minkowski term only, which may be called f AM ; 

fAM = _ e ±s 2 Ve. (8) 



By integrating this force across the boundary we obtain the vector surface force density 
" 2 f AM dn = ( x AM n, (9) 



n 



where a AM is the scalar 

1 



a 



AM 



-e (ei - e 2 J 



P2 , e J_c2 

e 2 



(10) 



when referring to the fields on the surface in medium 1. Here £t denotes the field 
component parallel to the surface, and 8^ the component normal to it. If reference is 
instead made to the fields in medium 2, 



a 



AM 



2 eol£l 



e 2 



E 2 + -E 2 

C-rp -|- Ojy 



The surface force, in general acting in the direction of the medium of lower permittivity, 
is thus in the present case directed along the normal vector n. It corresponds to <r AM 
defined as a positive quantity. 
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In accordance with usual conventions we let Sn denote the field component in the 
plane of incidence and £± the component normal to it. Thus in medium 1, 



El = £ 2 cos 2 9 t + £l, £ 2 N = £\\ sin 2 9 t , 
which together with Snell's law enables us to write the surface pressure as 

a AM = i( Cl - e 2 ) [(cos 2 9 t + sin 2 e t )£% + £l\. 



(12) 



(13) 



We now introduce the energy transmission coefficients Tji and T±, following the 
conventions of Stratton [28]. If £\ and £ l ± denote the components of the incident 
field in medium 1, and similar notations for the transmitted fields in medium 1, we have 

2 



Hi cos 9 t ( £\\ 

n 2 cos 9i \ £\ 

Til cos 9 t 

n 2 cos 9i \£ 




sin 29 i sin 29 1 



sm 2 (9 i + 9 t )cos 2 (9 i -9 t 



sin 26 i sin 26 t 

+ 9t)'- 



sm 



(14) 



(15) 



with Tii = a/^i, n 2 = y/t2. Let a denote the angle between £® and the plane of 
incidence, 



£u = £ l cosa, £\ 



£ l sin a. 



Then, by introducing the mean intensity I of the incoming beam, 



e n 2 c(£ l ), 



(16) 
(17) 



we have 



a 



AM I n l n 2 cos @i I-/ • 2 a i 2 /) \m 2 

[(sm 0j + cos Pt)Tj| cos a 



2c ni cos 6* t 

+ T ± sin 2 a]. (18) 

We assume circular polarization in the following, so that (sin 2 a) = (cos 2 a) = 1/2. 
It implies that also the hydrodynamical response of the surface becomes cylindrically 
symmetric. Using ordinary cylindrical co-ordinates we have dh/dd = 0, and sin#j = 
h r cosdi = h r (l + /i 2 ) -1 / 2 . Let now n denote the relative refractive index, 

n 2 



n 



< 1. 



m 



Then we can write 



a AM (r) 



2n 2 I(r) 1 — n 



f(n, h r ), 



c 1 + n 
where f(n, h r ) is the function 

1 + (2 - n 2 )h 2 r + hf + nhl^/l + {l-n 2 )h 2 r 



(19) 



(20) 



f(n, h r ) 



X 



nh 2 + VI + (1 - n 2 )h 2 
(1 + n) 2 



n 



+ y/1 + (1 - n 2 )h 2 



2 ■ 



(21) 
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We can now write the pressure condition at equilibrium as 



{p 2 -px)gh{r) - 

r dr 



rh r 



-a AM (r) 



(22) 



(recall that since the elevation is negative, h(r) < 0, whereas a (r) > when n\ > 712). 
The radiation pressure on the right hand side thus displaces the surface downward, while 
the surface tension term acts upward, opposing the depression. The sum of these two 
effects must be balanced by the influence of gravity (i.e., buoyancy) for mechanical 
equilibrium. Equation (1221) is our governing equation. Note that in this section the 
surface tension a has been assumed constant. 



3. First option: Investigation of a set of trial functions 

It is quite natural to check if the governing equation (|22|) possesses solutions 
corresponding to an abrupt jump from one kind of surface deflection to another, 
completely different one. If there are special solutions compatible with the protuberance 
form seen in the experiments it should be possible to retrace their form, at least 
approximately, by inserting reasonable test solutions into the equation itself. If we are 
on the right track, the difference between the right and left hand sides of the governing 
equation should be small. 

A number of trial functions were used, guided by experimentally observed shapes. 
Equations were put into non-dimensional form, introducing the Bond number Bq and 
the capillary length l c as 



2 



with u> the laser beam waist. Convenient non-dimensional (positive) variables for the 
height and radius of the deformation were 

H * fl(J0=rM (24) 

lc Wo 

Assuming a Gaussian laser beam, the intensity could be written in terms of the laser 
power, P, as 

J(r) = ^e~ 2R2 . (25) 

TTWq 

The governing equation could thus be expressed in non-dimensional form, with boundary 
conditions 

lim R = 00, lim Rh = —00. (26) 

We do not go into any detail concerning the use of these trial expressions; the main 
conclusion that we can make is that it is highly unlikely that there are special solutions 
of Eq. (122]) compatible with the observed form. This is consonant with the finding 
of Refs. [SSlEl], where a direct numerical simulation method is used. In no case did 
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the introduction of a "shoulder" shape improve the solution (i.e. tend to better equate 
the left and right hand side of Eq. (122]) ); always the contrary was observed. Although 
the result was negative, this brief discussion might be a useful inclusion to the future 
researcher of this question. 



4. Second option: Local variations of physical parameters. Heating effects 

Now focus attention on the physics of the situation: the central region of the laser beam 
where the field is strongest, may be expected to heat the fluids in the central region 
causing the physical parameters An, Ap and o to attain locally different values. A 
different shape can accordingly in principle be the result of the condition of mechanical 
equilibrium establishing a cylinder-like lower protuberance of radius ao and length h Q ; 
cf. Fig. [TJ The dimensions and other characteristics of the protuberance are determined 
by this equilibrium condition. 

In this section we shall use the values of ho and a$ measured in experiment, together 
with the known physical parameters at temperature 310. 65K and the temperature 
behaviour of An and o~o to estimate what local temperature near the tip might give 
rise to a protuberance of this smaller radius. The assumption made is thus that the 
heating by the laser is restricted to a small area near the tip of the deformation. We 
do not, therefore, include in our calculation the change to the density of the fluid in the 
boyancy term of the force balance, since that term is an integral over the full volume 
of the protuberance, of which only a small charge will be affected. The temperature 
change does, however, tend to increase the local surface tension cx(T) and the difference 
in index of refraction, An. The former increase tends to pull the protuberance upward 
(smaller value of h$) whereas the latter increases the electromagnetic radiation force 
and tends to pull downwards. 

Let us go through the details of the calculation. The mechanical forces acting on 
the protuberance region are as before from three different contributions: buoyancy (i.e. 
gravity or hydrostatic pressure), surface tension and radiation force. The pressure at its 
upper end is the hydrostatic pressure pigL, where L denotes the height of the depression 
at the laser power just before the protuberance becomes formed (cf. Fig. 1). Similarly, 
the pressure of the outside of the lower tip is also the hydrostatic pressure. At the lower 
end of the cylindrical section it is p2g{L + ho) (for simplicity we take all heights to be 
positive quantities). The buoyancy force on the protuberance is accordingly, when we 
include the volume of the hemisphere, 



B = Fb = g^o?o 



2 2 
(L + h + -a )Ap + pi(h + -a ) 



(27) 



where we have used p 2 = Pi + Ap. Subtracting off the weight W = gnalpi(h + |a ) of 
the liquid in the protuberance region we get the upward directed buoyancy force 

F B -W = gna 2 (L + h + ^a )Ap. (28) 
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Figure 3. Model geometry for estimation of tip temperature. 



Although the simplified test geometry thus implicitly assumed underestimates the 
volume of the deformation, we find that the results obtained are very insensitive to this. 
Increasing the volume by 50%, for example, only changes the estimated temperature 
increase by about 0.2K, which is not significant given the simplicity of the model itself. 

The buoyancy force has to be supplemented with the surface tension force 27ia a , 
also acting upwards. The third force, acting downwards, is the radiation force F ra d on 
the lower hemisphere tip. We shall work with the magnitude |-F ra d| of F rac j to avoid 
negative quantities. The condition for mechanical equilibrium of the protuberance can 
now be written as 

2 

gnal(L + h + -a )Ap + 2TTa a = |F rad |. (29) 



Once _F ra d is known - cf. the next section - Eq. (1291) determines the lower-tip surface 
tension o~q. Recall that the input parameters, to be inferred from experiments, are 
Ap, L, ho, and ao- 

We introduce now a quantity Q, as a non-dimensional measure of the force exerted 
by the incident laser beam. It is defined as follows: If the beam were plane, with 
intensity Jo, we would according to common usage write the radiation force F ra ^ (here 
taken positive) on a cross-sectional area *na\ as 

-Frad = ^^-nalQ, (30) 

meaning that a perfect absorber corresponds to Q = 1. Considering instead a Gaussian 
beam with total power P and beam waist we have, if Iq now means the intensity on 
the symmetry axis, 

T 2n 2 P 



TTWq 



(31) 
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Accordingly, we can in the present case write 



_ InlPal 

rad 



2 °Q. (32) 



rad 



^.1. Electromagnetic force on a hemisphere 

As is known, the force acting on a closed surface S may be calculated by integrating 
the stress tensor 

T = S®T> + n®B-\{S-T> + n-B)l (33) 
over the surface 

<J) dSr ■ (T) ■ z = I dSr ■ (T rz ) ■ z (34) 
Js Js 

with 1 the unit matrix and (• • •) denotes time average and hats denote unit vectors. 
(For an isotropic, dielectric medium the tensors of Minkowski and Abraham coincide. 
The Abraham force [9] oscillates out and gives no contribution in the optical case.) Here 

£ = Re{Ee iaJ '}, etc., (35) 

where E, D,B and H are complex field vectors. For two field components X i7 Xj we 
have (XiXj) = \Re{X t X*}. 

In the case of a full sphere, or indeed the total force on any isolated body in 
a homogeneous external medium, it may be opportune to integrate the stress tensor 
over a closed surface far from the body in order to make use of simplified asymptotic 
expressions for the spherical Bessel functions involved. When evaluating the force on a 
part of the body surface, however, integration must be performed at the actual interface, 
and we find it simpler in this case to express the stress tensor both outside and inside 
the spherical surface in terms of the interior fields, since the external fields have both 
incident a scattered components. Using the continuity of D t ^Eq^E^ and H across the 
surface, the axial front force on a sphere segment # o <$< 7r ) O<0<27r can be 
expressed as 

Frad = 4 dcj) sin 9d6{T r e ? - Tf?) 
Jo Je 

2 2 fir 
7m 2 a £ /-2 t\ / /-2|7TW|2 , lcnu|2\ 



n 2 -l) / (n 2 |£ r T + |£ t T) 







2 

xcos0sin0d0 (36) 

where \E t \ 2 = l-E^I 2 + l^l 2 ? an d we have used the fact that the integrand has no 
dependence. Superscript w signifies that the fields are evaluated just within the surface 
of the hemisphere, at radius r = a~~ . The expression is obviously real and there is no 
longer a need to explicitly take the real part. 

We now make the assumption that the electromagnetic fields inside the 
hemispherical surface equal those inside a full sphere, and use the electromagnetic field 
expressions due to Barton and co-workers (31], which we quote in some detail in the 
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Appendix. As before, we let the incident field be a circularly polarized plane wave, 
for which the explicit field expansions are found in Eqs. (IA.12j) - flA.14j) . and we assume 
the laser width wq to be sufficiently large that laser intensity can be approximated as 
uniform over the hemisphere (this is a reasonable approximation for our purposes since 
the optical force on a hemisphere is quite concentrated near the symmetry axis [13J .) . 
We make the convenient definitions 



i' +2 [n^™)£P(«) " Vftna)^]" 1 (37) 
d l = i< +1 [V^a)£p(a) - ^'Neft 1 (38) 
where ipi and £/ are Riccati-Bessel functions of order /, 

a = kr= 2 ^l = ^l (39 ) 
A 2 c 

is the number of hemisphere circumferences per optical wavelength (in medium 2). Thus 
we find the radiation force on the sphere segment to be 

— 2 oo oo 

Q = n ^ 4 ^2^2^c k c1ij k (na)iJi(na)I k i 

k=l 1=1 

+ o 2 [cjkC*'0fc(^ Q! )'0[(^a) + d k d^ip k (na)ipi(na)}M k i 

+ a 2 [c k d^ k (na)ijji(na) + d k c*^ k (na)^(na)]N k i \ (40) 



where the coefficients I k i, M k i and N k i contain the actual polar angle integration over 6 
from 7r/2 to 7r and are given as 



iki 



(2k + 1)(2Z + 1) J° duuP k \u)Pl(u), (41) 



2k + 1 21 + 1 f° duu fnU . 



k(k + 1)1(1 + 1) J_ x 1-u 2 
+ (l-u 2 fPl'(u)Pl'(u)l (42) 

^-- ^g^ /^^- (43) 
Numerically, these coefficients may be tabulated once and for all, and calculation is thus 
not particularly expensive. One finds that the sums in ( HOj) can be truncated at a value 
somewhat greater than a. In our case a is in the order of 60, so about 4000 terms were 
evaluated to calculate Q. 



4-2. Numerical results 

Given the value of Q calculated from f PfOl) using the unperturbed An, using the fact that 
Q oc An(T) and inserting the temperature dependences for An, Ap and cro from Eqs. (|2J)- 
(jH), temperature T is the only unknown in the equation of mechanical equilibrium, and 
may thus be determined from Eq. fl29|) . Explicitly we write 

B(T ) + S(T )e 2 » - |F rad |(T )6 /3 = (44) 
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with B the buoyancy force, S the contribution from surface tension [see eq. (j29p ]. and 
6 = ^^; T = Tc + e(T -Tc) (45) 

with T = 310.65, the external temperature. The equation is easily solved with respect 
to O using Newton's method. 

Table [1] shows values of <To and An, and corresponding tip temperatures T new , in 
the cases for which experiments were performed OUSIHH] • There are several conclusions 
to be drawn from these values: 

(i) It is apparent that cr is increased in comparison with the initial global value 
o = 2.413 x 1(T 7 N/m, calculated from Eq. (0}. 

(ii) Correspondingly, the new tip temperature T new is also increased. There occurs 
thus a heating of the tip, in accordance with the expectation. It should here be 
noted that the relationship between surface tension and temperature is in our case 
counterintuitive: increasing surface tension means increasing temperature, in the 
region above the critical point. 

(hi) The increase in temperature is relatively large; for moderate powers P the values 
of T new lie about 4 degrees higher than the ambient temperature T . It may be of 
interest to compare this with the much smaller temperature increase occuring in a 
homogeneous fluid illuminated by a laser beam. Adopting the Gaussian form for 
I(r) from Eq. (|25|) we may solve the heat conduction equation 

V 2 T(r) + -J(r) = (46) 

K 

in cylindrical symmetry, where a ~ 3 x 10 -4 cm^ 1 is the thermal absorption 
coefficient and k = 1.28 x 10~ 3 W cm -1 K _1 is the thermal conductivity. One 
finds [31 [23] that on the symmetry axis r = the local temperature increase is 
AT « aP/(47T«)(ln 100 + 7), where 7 « 0.5772 is Euler's constant. With P = 1 W 
this amounts to about 0.1 K, which is a very moderate heating. 

The difference in expected temperature increase mentioned in point (lm|) is not 
as unreasonable as it seems. Our calculation concerns a local effect, the heating of 
the interface near the protuberance tip and its immediate surroundings, whereas the 
calculation in point flmj) is a heating of a much larger volume of fluid. To heat a finite 
body of fluid by the same extent, would require a much higher power. A significantly 
larger heating is expected near the tip of the protuberance anyway, since the deformation 
acts as a lens, focussing the incoming light and giving rise to local intensity maxima 
much higher than that near a laser beam in a homogeneous medium (see e.g. [30]). 

To further illustrate the importance of local geometry, consider the following simple 
argument: Let the interface be modeled as a thin horizontal plate of thickness Az and 
surface area A = no 2 , surrounded by a vacuum, illuminated by a laser power P. When 
thermal equilibrium is established, the rate of absorbed heat aP has to balance the 
rate of radiated energy to both sides, i.e. 2Ao~sb(T 4: — TJ 4 ), where o~sb is the Stefan- 
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w [fim] 


P [mW] 


L + h \pm\ 


a [p,m] 


Q [io- 3 ] 


An 


ao [10- ? N/m] 


T new [K] 


6.3 


1200 


70 


3.2 


4.564 


0.0109 


15.7 


319.1 


6.3 


600 


40 


3.8 


4.562 


0.0090 


7.57 


314.3 


4.8 


1200 


72 


2.5 


4.567 


0.0121 


24.0 


323.6 


4.8 


600 


43 


2.6 


4.567 


0.0096 


9.80 


315.7 



Table 1. Parameter sets from experiments [3j[16j[T8] and calculated value of non- 
dimensional radiation force Q from these. The requirement of mechanical equilibrium 
allows us to estimate the new tip temperature T ncw and locally diffcnt values of An 
and Co- Parameter values are g = 9.81m/s 2 , Ap = 59.61kg/m 3 , To = 310. 65K, u(To) = 
2.413 x lCr 7 N/m, n 2 = 1.46, An(T ) = 0.00672. 



Boltzmann constant. Since the excess temperature AT = T — To is small, AT /To <C 1, 
we get the relationship 

AT = aP ^ . (47) 

Inserting a as above, and with P = 1 W, <t S b = 5.67 x 10~ 8 Wm" 2 K~ 4 , T = 300 K, we 
get 

AT = (7.80 x 10~ 4 Km)^. (48) 

cr 

Now take for defmiteness a = 4 /im. If we choose the thickness to be small, Az = 1 nm, 
we get only a small temperature increase, AT « 0.05 K. With Az = 100 nm the result 
is much higher, AT w 5 K. We thus see that the thermal behaviour is highly sensitive 
to the thickness over which absorption takes place. In other words, local temperature 
variations can take very different values from from overall (global) ones. 



5. Concluding remarks 

We thus arrive at an explanation of the observed effect that seems quite plausible, 
namely that the tip is heated locally and so makes the interface thermodynamically 
nonuniform. The change in physical parameters An and o~o conspire to favour a narrower 
deformation. We have performed a fairly simple model calculation which indicates that 
this is the case, yet does not constitute a full explanation for the particular shape that 
the deformation takes. 

It seems likely, however, that the increase in surface tension in the areas where 
the most laser light is absorbed, i.e. at the tip of the emulsion and at the "shoulders" , 
could relate directly to the onset of a Plateau-Rayleigh-like instability [33]. A Plateau- 
Rayleigh instability is driven by surface tension trying to decrease the interface area, 
and surface tension is increased whereever the liquids are heated. As we have discussed 
above, a deformed surface will act as a lens, causing local light intensities to far exceed 
average ones near areas of sharp deformations, such as "shoulders" and tip. A slight 
change in shape increases the heating due to light focussing, increasing the surface 
tension and giving further change in shape, and so on. In this way, local heating can 
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be thought to drive the required instability. Instead of a separate droplet forming (as 
indeed it does for even higher laser powers [15]), the radiation pressure may serve the 
role of stabilizing the intermediate "shoulder" geometry, similar to the situation reported 
in [17] . A further investigation into this stability issue is certainly warranted, and of 
great potential interest. 

The other option considered in this paper, namely that the shape can be explained 
merely by balancing the forces of laser pull, buoyancy and surface tension with 
unchanged physical parameters, turned out not to be supported. Our extensive 
numerical search in this direction led to no indication that the unperturbed force balance 
equation has solutions at all similar to the observed deformations. Of course, trials of 
this sort cannot lead to a decisive falsification. Nevertheless, we feel our investigation 
supports a firm conviction that that the observed deformations have their roots in local 
effects connected with temperature variations, and cannot be due to global effects. 

Finally, we suggest that the proposed explanation should be possible to confirm 
experimentally, were the experiment [15] to be set up again. We envisage that use of 
an infrared (thermal) camera could enable detection of local temperature gradients. A 
more complete theory could moreover be assisted by a simulation in which absorptive 
heating were taken into account, although this would be a project suited for heavy 
numerics rather than analytical means. In light of the large and growing interest in 
fluid manipulation with laser light, however, understanding the instability studied herein 
could be of some importance. 

Appendix A. General formulae 

Appendix A.l. From Mie theory 



The internal electric field components are expanded in spherical harmonics according 
to [31] 



oo 





1=1 m=—l 



oo 




1=1 m=—l 




(A.2) 



1=1 m=—l 




(A.3) 



where Vt = 9,4> and the coefficients 




(A.4) 
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and the incident field E* is contained in the quantities 
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(A.5) 

B lm l(l + l)^(a)J Hi/ H (A ' 6) 

where incident fields are evaluated at r = a + and the integral is over all solid angles. 
Here H = £ /(c/i ). 

Appendix A. 2. Circularly polarised plane wave 

A circularly polarised plane wave propagating along the z direction may be expressed 
as ( [32], section 10.3) 



E* = E Ae lkz = E J2 WMU + WiMXn 



i=i 



+ _v x ji{kr)X n ] 



with 



Xj[ m (n) 



i 



where k = n 2 uj/c and A = x + iy. The radial component may then be written 
section 10.4) at r = a + 

E t = Hf E i ,+ V4ttW + 1)(2Z + 



(A.7) 
(A.8) 



(A.9) 



The radial magnetic component is now found from Maxwell's equations as 

in 2 



if! 



// c 



By using formulas (IA.6[) and (I A . 9 1) we find 

/4tt(2/ + 1 



.4 



■-> \ I j i j I, "ml -^im 
" > /(/ + 1) 11-2 



and the field components (1A. 1|) . (1A.2I) and (1A.3[) may be written 
= ^(2/ + I)q^(^)^ 1 
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(A.10) 
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(A.13) 



+ ditjji(na)Pf sin6 



sin# 
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having suppressed the argument cos# of the Legendre polynomials. Coefficients q and 
d\ are defined in Eqs. fl37|) and fl38|) . 
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